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Abstract Because the propagation of neutrinos is affected by the presence 
of Earth matter, it opens new possibilities to probe the Earth's interior. Dif- 
ferent approaches range from techniques based upon the interaction of high 
energy (above TeV) neutrinos with Earth matter, to methods using the MSW 
effect on the oscillations of low energy (MeV to GeV) neutrinos. In princi- 
ple, neutrinos from many different sources (sun, atmosphere, supernovae, 
beams etc.) can be used. In this talk, we summarize and compare different 
approaches with an emphasis on more recent developments. In addition, we 
point out other geophysical aspects relevant for neutrino oscillations. 

Keywords Neutrino absorption • Neutrino attenuation • Neutrino oscilla- 
tions • Matter effects • Neutrino tomography 



Supported by the W. M. Keck Foundation and NSF grant PHY-0503584. 
Walter Winter 

Institute for Advanced Study 
School of Natural Sciences 
Einstein Drive 
Princeton, NJ 08540 
E-mail: winter@ias.edu 



2 



Walter Winter 



1 Introduction 

Neutrinos are elementary particles coming in three active (i.e., weakly in- 
teracting) flavors. Since the cross sections for neutrino interactions are very 
small, neutrinos practically penetrate everything. However, one can com- 
pensate for these tiny cross sections by just using enough material in the 
detector. Depending on neutrino energy and source, the detector has to be 
protected from backgrounds such that the neutrino events cannot be easily 
mixed up with different particle interactions. Neutrinos are produced in de- 
tectable numbers and with detectable energies by nuclear reactions in the 
sun, by cosmic ray interactions in the Earth's atmosphere, in nuclear fission 
reactors, in supernova explosions, in the Earth's crust, and possibly by astro- 
physical sources. In addition, accelerator-based neutrino sources specifically 
designed to produce a high-intensity neutrino beam have been successfully 
operated (such as K2K or MINOS) or are planned. Thus, there are neu- 
trinos from various different sources with different energies. 

One of the most recent exciting discoveries in neutrino physics is neu- 
trino oscillations, i.e., neutrinos change flavor while traveling from source to 
detector. This quantum mechanical phenomenon implies that neutrinos mix, 
i.e., the eigenstates of the weak interaction are not the same as the mass 
eigenstates, and at least two out of the three have non-vanishing masses. 
This is probably the most direct evidence today for physics beyond the stan- 
dard model of elementary jjarticle physics. Recent neutrino oscillation ex- 
periments, especially SNO l], KamLAND jl^l, Super-Kamiokande 24], and 
CHOOZ j2] have helped to quantify this picture. Unlike in quark mixing, two 
out of the three mixing angles are large, and one is even close to maximal. In 
addition, the oscillation frequencies have been fairly precisely measured. For 
one of the mixing angles ^13, however, only an upper bound exists, and several 
parameters (the arrangement of masses, i.e., mass hierarchy, and one com- 
plex phase Sep relevant for neutrino oscillations) are still unknown. Future 
experiments will probe these parameters starting with the Double Chooz , 
T2K ^3 , and NO^A [ll| experiments (for the prospects for the next decade, 
see, e.g., Ref. ^^). 

For neutrino tomography the relevant aspect is the sensitivity to Earth 
matter. Since it is well known that the cross sections with matter rise at least 
until lOTeV ^6Q] . the probability of matter interactions can be increased by 
higher neutrino energies. Neutrino absorption tomography uses this effect to 
infer on the matter structure. For neutrino oscillations, we know that the so- 
called MSW matter effect 48; 49] is the most plausible explanation for the 
solar neutrino deficit |2lj ]. This, however, implies that neutrino oscillations in 
the Earth have to experience this effect, too. Neutrino oscillation tomography 
uses the MSW effect to study the matter structure. 

2 Tomography using the propagation of neutrinos 

Tomography using the propagation of neutrinos [sit l50| assumes a neutrino 
source with a well-known flux and flavor composition, a well-understood neu- 
trino detector, and a speciflc neutrino propagation model between source and 
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detector. The key ingredient to any such tomography is a considerable depen- 
dence of the propagation model on the matter structure between source and 
detector. Compared to the detection of geoneutrinos, the object of interest is 
not the neutrino source, but the material along the baseline (path between 
source and detector). If the matter structure along the baseline is (partly) 
unknown, the information from counting neutrino events at different ener- 
gies by the detector can be used to infer on the matter profile. Two accepted 
propagation models could be used for neutrino tomography: 

Neutrino absorption: Because the cross section for neutrino interactions in- 
creases proportional to the energy, neutrino interactions lead to atten- 
uation effects. Useful neutrino energies for a significant attenuation are 
> 1 TeV. 

Neutrino oscillations: The MSW effect in neutrino oscillations 

(coherent forward scattering in matter) leads to a relative phase shift of 
the electron flavor compared to the muon and tau flavors. This phase 
shift depends on the electron density. Useful neutrino energies require 
substantial contributions from the MSW effect as well as large enough 
oscillation amplitudes. Depending the relevant Am?, neutrino energies 
between 100 MeV and 35 GeV are optimal for studying the Earth's inte- 
rior. 

Beyond these two models, at least small ad-mixtures of non-standard effects 
have not yet been excluded. Some of these non-standard effects are sensi- 
tive to the matter density, too. Examples are mass- varying neutrinos with 
acceleron couplings to matter fields [43, non-standard neutrino interactions 
(see Ref. jsJl and references therein), and matter-induced (fast) neutrino de- 
cay 28]. Because there is not yet any evidence for such effects, we do not 
include them in this discussion. 

Given the above neutrino energies, there are a number of potential sources 
which could be used for neutrino propagation tomography. For neutrino oscil- 
lations, solar neutrinos, supernova neutrinos, atmospheric neutrinos, and neu- 
trino beams (such as superbeams or neutrino factories) are potential sources. 
For neutrino absorption, high-energy atmospheric neutrinos, a possible high- 
energy neutrino beam, or cosmic sources are possible sources. 

As far as potential geophysics applications are concerned, neutrinos may 
be interesting for several reasons: 

1. Neutrinos propagate on straight lines. The uncertainty in their path (di- 
rection) is only as big as the surface area of the detector. 

2. Neutrinos are sensitive to complementary quantities to geophysics: Neu- 
trino absorption is directly sensitive to the matter density via the nucleon 
density. Neutrino oscillations are sensitive to the electron density which 
can be converted in the matter density by the number of electrons per 
nucleon (for stable "heavy" materials about two). On the other hand, 
seismic wave geophysics needs to reconstruct the matter density by the 
equation of state from the propagation velocity profile. 

3. Neutrinos are, in principle, sensitive to the density averaged over the 
baseline, whereas other geophysics techniques are, in principle, less sensi- 
tive towards the innermost parts of the Earth. For example, seismic shear 
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waves cannot propagate within the outer hquid core, which means that a 
substantial fraction of the energy deposited in seismic waves is reflected 
at the mantle-core boundary. Other direct density measurements by the 
Earth's mass or rotational inertia are less sensitive towards the innermost 
parts, too, because they measure volume-averaged quantities. 

Given these observations, there may be interesting geophysics applications 
exactly where complementary information is needed. Possible applications 
range from the detection of density contrasts in the Earth's upper or lower 
mantle, to the measurement of the average densities of the outer and inner 
core by independent methods. 



3 Neutrino absorption tomography 

Here we discuss tomography based on attenuation effects in a neutrino flux 
of high enough energies, which we call, for simplicity, "neutrino absorption 
tomography" . After we have introduced the principles, we will discuss pos- 
sible applications with respect to tomography of the whole Earth as well as 
specific sites. 



3.1 Principles 

"Neutrino absorption tomography" uses the attenuation of a high-energy 
neutrino flux as a propagation model. In this case, weak interactions damp the 
initial flux by the integrated effect of absorption, deflection, and regeneration. 
For example, muons produced by a muon neutrino interaction are absorbed 
very quickly in Earth matter, whereas tauons produced by tau neutrinos 
tend to decay before absorption (and some of the decay products are again 
neutrinos). Only the integrated effect leads to attenuation of the flux. The 
magnitude of the attenuation effect can be estimated from the cross section 

^ ^ 10-35 (1) 

E TeV ^ ' 

to be of the order of several per cent over the Earth's diameter for E„ = 
ITeV. The interaction cross section rises linearly up to about 10 TeV |60|. 
whereas the behavior above these energies is somewhat more speculative. The 
energies are usually as high as standard neutrino oscillations do not develop 
within the Earth. Since the neutrinos interact with nucleons, the attenuation 
is directly proportional to the nucleon density. Therefore, neutrino absorp- 
tion is a very directly handle on the matter density with an extremely tiny 
remaining uncertainty from composition and the difference between neutron 
and proton mass. 

As far as possible potential neutrino sources are concerned, Eq. |^ re- 
quires very high neutrino energies. The existence of corresponding neutrino 
sources is plausible and will be tested by upcoming experiments commonly 
referred to as "neutrino telescopes". These neutrino telescopes could also 
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a) Isotropic flux 



b) High-energy 
neutrino beam 



c) Cosmic point 
source 



Fig. 1 Three different approaches to "Whole Earth Tomography" using neutrino 
absorption. The hues refer to different basehnes. 

serve as prototypes for the detectors useful for neutrino absorption tomog- 
raphy. The only detected source so far is atmospheric neutrinos produced by 
the interaction of cosmic rays in the atmosphere. Unfortunately, the atmo- 
spheric neutrino flux drops rapidly with energy, which means that statistics 
is limited at the relevant energies > 1 TeV (see, e.g., Ref. for specific 
values). Other potential candidates include many different possible astro- 
physics objects, as well as particle physics mechanisms such as decays of 
dark matter particles. Since we know that the Earth is hit by cosmic rays of 
very high energies, it might be inferred that astrophysical mechanisms exist 
which accelerate particles (for example, protons) to these high energies. It 
is plausible that such mechanisms also produce neutrinos. Potential mecha- 
nisms could either produce discrete fluxes from individual objects, or their 
integrated effect could lead to a diffuse flux over the whole sky. Eventually, 
one could think about a neutrino beam producing high-energy neutrinos. If, 
for instance, one used the protons from LHC (7 TeV) to hit a target, the 
decaying secondaries (pions, kaons) would produce a neutrino flux peaking 
at about 1 TeV. 



3.2 Whole Earth tomography 

For possible applications of neutrino absorption tomography, there exist two 
different directions in the literature: Either one could "X-ray" the whole 
Earth ("Whole Earth tomography"), or one could think about the investiga- 
tion of specific sites in the Earth's mantle. We summarize in Fig. ^ different 
approaches to "Whole Earth tomography" . In case a) (isotropic flux) a neu- 
trino flux from many directions is detected by a detector with good directional 
resolution. For instance, a possible neutrino source would be a cosmic diffuse 
flux ^3 (related work: Ref. |6j ) or the high-energy tail of atmospheric neu- 
trinos (see, e.g., Ref. H^). This application could be very interesting because 
it might be available at no additional experimental effort. However, if one 
wants to study the innermost parts of the Earth, it is (except from sufficient 
directional resolution and flux isotropy) a major challenge that the fraction 
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Fig. 2 Different possibilities for neutrino tomography using a high-energy neutrino 
beam. The labels "IP" refer to possible interaction points. See text for more details. 



of the sky which is seen through the Earth's inner core is very small (~ 1%), 
which means that the statistics for this specific goal is very low. Very good 
precisions may, on the other hand, be obtained for the mantle (see Fig. 4 
of Ref. ^3)- (high-energy neutrino beam) the 

detector is moved to obtain many basehnes, whereas the source is kept fixed. 
In this case, high precisions could be obtained However, a major chal- 
lenge might actually be the operation of a high-energy neutrino beam with 
a moving decay tunnel. Note that such a beam could not only be used for 
whole Earth tomography but also for local searches (see below). In case c) 
(cosmic point source) [M: l45|. the flux from a single object is used for the 
tomography of the Earth. In this case, the flux has to be constant in time to 
be detected either by a moving detector, or by one detector using many base- 
lines by the rotation of the Earth. Note that the second mechanism cannot 
be used for the currently largest planned neutrino telescope "IceCube" 
because it is residing at the south pole. 



3.3 Specific site tomography 

Compared to "whole Earth tomography" , a different direction is the investi- 
gation of individual sites, such as in the Earth's mantle. For example, Ref. 
extensively reviews techniques based on a high-energy neutrino beam. We 
summarize some of those in Fig. |21 The neutrinos, produced by the source 
"S" , may interact at several possible interaction points IP. If, for example, 
the site of interest is the dark-shaded cavity, an interaction at IPi could 
create a particle shower leading to sound production, which may be detected 
by a microphone array at the surface. In addition, the final neutrino flux 
detected at "D" would be damped depending on the material density in the 
cavity. An interaction at IP2 just below the surface (< 200m) would pro- 
duce muons which could still be detected at the surface (such as possibly by 
a muon detector on a truck). A variation of this flux detected by a moving 
muon detector could point towards heavy materials. Eventually, a neutrino 
interaction at IP3 within the sea water below a muon moving detector would 
indicate that the initial neutrino has arrived. Since the neutrino energy de- 
creases rapidly by moving the detector out of the beam axis by kinematics, 
attenuation effects also decrease and the initial flux could be measured by 
the "off-axis" technology. Comparing this flux to the on-axis flux reveals the 
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attenuation along the path and therefore some information on the matter 
structure. 

In summary, there are many potential applications of neutrino absorp- 
tion tomography. The coming years, especialy the operation of IceCube, will 
reveal the possible existence of high-energy cosmic neutrino fluxes. Operat- 
ing a high-energy neutrino beam may be a major technical challenge, which 
definitively needs further investigation. 

4 Neutrino oscillation tomography 

In this section, we discuss neutrino tomography using oscillations. First, we 
introduce the principles of neutrino oscillation tomography: Neutrino oscil- 
lations in vacuum and matter, numerical approaches to neutrino oscillation 
tomography, as well as conceptual (mathematical) problems. Then, we show 
applications related to solar and supernova neutrinos, and we discuss tomog- 
raphy with neutrino beams. 

4.1 Principles 

Neutrino oscillation tomography uses neutrino oscillations in matter as prop- 
agation model. Possible neutrino sources include "natural" ones {e.g., sun, 
supernova, atmosphere), as well as "man-made" ones {e.g., superbeam, /3- 
beam, neutrino factory). The detection technology depends on the neutrino 
energy and ranges from water Cherenkov detectors (lower energies), over liq- 
uid scintillators (medium energies), to iron calorimeters (high energies), just 
to mention some examples. 

Neutrino oscillation phenomenon 

Neutrino oscillations are a quantum mechanical phenomenon with two pre- 
requisites: First, the weak interaction eigenstates have to be different from the 
propagation/mass eigenstates (flavor mixing). Second, the neutrino masses 
have to be different from each other, which implies that at least two of the 
active neutrinos have to have non-zero mass 12]. In the limit of two flavors, 
the flavor transition probability Va — > i^p in vacuum can be written as 



where 9 is the mixing angle of a 2 x 2 rotation matrix U, Am? = — 
is the mass-squared difference describing the oscillation frequency, L is 
the baseline (distance source-detector), and E is the neutrino energy. Note 
that the quotient L/E determines the oscillation phase. Similarly, the flavor 
conservation probability z^Q. Va is given by Pcq = l—Pap from conservation 
of unitarity. Practically, Pap is measured as function of E (convoluted with 
the neutrino flux and cross sections) for a fixed baseline since the detector 
cannot be moved. Since we do know that we deal with three active flavors. 
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the complete picture is somewhat more complicated. Three-flavor neutrino 
oscillations can be described by six parameters (three mixing angles, one 
complex phase, and two mass squared differences), which decouple into two- 
flavor oscillations, described by two parameters each, in certain limits (see 
Ref. for a recent review). In summary, we have two almost decoupled 
two-flavor oscillations described by two very different frequencies and large 
mixing angles, often referred two as "solar" {Am^i, O12) and "atmospheric" 
(Z\m|j^, ^23) oscillations. Those could be coupled by ^13, for which so far 
only an upper bound sin^(20i3) < 0.1 7] exists. In addition, we do not 
yet know anything about the complex phase Sqp^ which could lead to sub- 
leading effects, and the sign of Am^i ("mass hierarchy"). These parameters 
will be probed by neutrino oscillation experiments in the coming years. In 
this section, we concentrate on the two-flavor case for pedagogical reasons. 



Matter effects in neutrino oscillations 



Key in gred ient to neutrino tomography are matter effects in neutrino oscil- 
lations [Ggliali^l. Since Earth matter contains plenty of electrons, but no 
muons or tauons, charged-current interactions of the electron neutrino flavor 
through coherent forward scattering lead to a relative phase shift compared 
to the muon and tau neutrino flavors. In the Hamiltonian in two flavors, the 
matter term enters as the second term in 

\.^fA{n^)0 



in flavor space, where vl(rie) — -ii\/2G prie is the matter potential as function 
of the electron density Ue and the coupling constant Gf, and the different 
signs refer to neutrinos (plus) and antineutrinos (minus). Assuming that the 
number of electrons per nucleon is approximately 0.5 for stable "heavy" (con- 
siderably heavier than hydrogen) materials, the electron density can be con- 
verted into the matter density as — 0.5 p/mj^ with m^r the nucleon mass. 
In this case, there is some material dependence of this factor 0.5 ("electron 
fraction"), which, however, might also be used to obtain additional infor- 
mation on the composition. In two flavors and for constant matter density, 
Eq. 121 can be easily re-written by a parameter mapping between vacuum 
and matter parameters: 

Pa0=sin2(20)sinM^^j , (4) 



where 



with 



/Im^ = e • Am^ , sin(20) = ^^^^ , (5) 



^ EE ^sin^(2e') + (cos(2e') - Af , (6) 
^^2V2G^_ (7) 
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Fig. 3 Example for a REM- ("Reference Earth Model") based matter density 
profile for a baseline of 12 000 km as used for the numerical evaluation of the neu- 
trino propagation (20 steps). The matter density within each layer is assumed to 
be constant. 

One can easily read off these formulas that for A cos(20) the parameter 
^ in Eq. © becomes minimal, which means that the oscillation frequency 
in matter becomes minimal and the effective mixing maximal (c/., Eq. l(SJ)). 
This case is often referred to as "matter resonance", where the condition 
A cos(20) evaluates to 

Am'^ feV^l 

Src. ~ 13 200 cos(20) 7/ 31 • (8) 
p [g/cm'^J 

This condition together with the requirement of a large oscillation phase 
sin^(Z\m^L/(4i?)) = C(l) leads to the "ideal" energies for neutrino oscilla- 
tion tomography depending on the considered Am?: 

Aml^: lOOMeVtolGeV, 
Aml^ : E ^few GeV to 35 GeV . 

If the neutrino energy is far out of this range, either the matter effects or 
the overall event rate from oscillations will be strongly suppressed. However, 
there are also possible applications. Since, for instance, for solar neutrinos 
E E,-cs , one can use the absence of the resonance for analytical simplifica- 
tions, as we will discuss later. 

Numerical evaluation and conceptual problems 

In order to numerically study neutrino oscillation tomography, a commonly 
used method is the "evolution operator method" (c/., e.g., Ref. This 
method assumes that the matter density profile be discretized into layers 
with constant density (c/.. Fig. |21for an example). The initial state is 
then propagated through the different matter density layers with depths xj 
with the evolution operators 



Vixj,pj) = e-''^^P^^'=^ 



(9) 
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Fig. 4 Neutrino oscillation tomography using solar neutrinos for the investigation 
of a cavity in the Earth's mantle. 

while the Hamiltonian within each layer Ti. {cf., Eq. J^Jl) is assumed to have no 
explicit time-independence (it is given in constant density pj). The transition 
probabihty is then obtained as 

V(a;l,pl)|^.„)|^ (10) 

In practice, Eq. (|10|) is evaluated by diagonahzing the Hamiltonian for each 
density step, i.e., by calculating the mass eigenstates in each matter layer. 
Note that in general 

[Vix„p,),Vixj,pj)]^0 for p.^pj, (11) 

which means that the evolution operators of different layers do not necessarily 
commute. This already implies that the information from a single baseline 
must be somehow sensitive towards the arrangement of the matter density 
layers. This is very different from X-ray or neutrino absorption tomography 
which do not have positional information from one baseline. 

An important conceptual problem in neutrino oscillation tomography is 
the matter profile inversion problem [Tol : l20l | . Assume that a matter density 
profile such as in Fig. |31is given. For a specific experiment setup, it is then 
fairly easy to compute the corresponding transition probabilities or event 
rates as function of energy. However, the reverse problem is theoretically 
generally unsolved: Assume that the transition probability is known up to 
infinite energies, then it would be very useful to be able to compute the 
matter profile from that. So far, there have been several attempts to solve 
this problem using simplifications, such as 

— Simple models using only very few discrete steps (see, e.g., Refs. |E2IE3 

HI) _ , 

— Linearization in a low density medium (solar, supernova neutrinos) |j] 

— Discretization of a more complex profile using non-deterministic algo- 
rithms to fit a large number of parameters [ssf . 

Below, we will discuss some of these approaches in greater detail. 

4.2 Neutrino oscillation tomography with solar and supernova neutrinos 

Solar and supernova neutrinos are theoretically very interesting for neutrino 
tomography because matter effects are introduced off the resonance in Earth 
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matter, i.e., the neutrino energy E <C -Brcs (c/., Eq. for Am^ = Arn^i), 
or equivalently A ^ I. This means that one does not expect strong matter 
effects in Earth matter as opposed to within the sun. However, this hmit 
is theoretically very useful to study tomography because it allows for per- 
turbation theory and other simplified approaches. It is often referred to as 
neutrino oscillations in a "low density medium" [s^ |3^ because the density 
in the Earth is much lower than in the sun. 



Detecting a cavity 

We show in Fig.0]a possible setup for neutrino tomography using solar neu- 
trinos following Ref. [37'|. In this setup, the detector is fixed while the Earth is 
rotating, which means that the cavity with density p is "exposed" (in line of 
sight sun-detector) a time < toxp < 24 hr per day. The change in the oscilla- 
tion probability during this time is, depending on geometry and density con- 
trast, < 0.1%. This leads to a required detector mass M > 130 Mt/(icxp[hr]), 
which has a lower limit of 5Mt at the poles. Thus, from the statistics point 
of view, this approach is very challenging, and backgrounds might be an im- 
portant issue. In addition, for such large detectors, the detector surface area 
might be of the order of the cavity size. There are, however, interesting the- 
oretical results from such a discussion. Let us define the oscillation phases in 
the individual steps xi, d, and X2 as |37j 

^ ^ll^ 7sin2 20i2 + (cos20-i.)2 (12) 

with the corresponding matter potentials Ai {cf., Fig.^J). One can show that 
if mass eigenstates arrive at the surface of the Earth (solar and supernova 
neutrinos), the change in probability AP (cavity exposed-not exposed) only 
depends on ^2, but not on In addition, there is a damping of contributions 
from remote distances X2 , which means that solar neutrinos are less sensitive 
to the deep interior of the Earth than to structures close to the detector. 



Matter density inversion problem 

A further application of the low density limit is to theoretically solve the 
matter profile inversion problem. Following Ref. 0| , the Earth matter effect 
on solar or supernova neutrinos is fully encoded in the quantity ( "day- night 
regeneration effect" ) 



pnrght _ pday ^ _ ^^^2 ^^^2 20^^^(J) (13) 
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with 



f{S) = J dxA{x) sin 




2 / uj{x')dx' 



J {S cos 26*12 - A{x) /2)2 + sin^ 26*12 : 



A{x) ^ V2GFne{x) , 5 : 



AE 



(14) 

(15) 
(16) 



This imphes that the measured quantity is i.e., a function of energy, 

which needs to be inverted into the matter profile A{x). EspeciaUy, the double 
integral in Eq. H14(l is quite complicated to invert. However, using the low 
density limit A <^ 25 (or equivalently A <C 1) as well as AL <C 1 (i <C 
1 700 km) , one can linearize Eq. H14|) in order to obtain 



f{5) = J dxA{x) sm[2S{L - x)] . 



This is just the Fourier transform of the matter density profile, i.e., 



(17) 



A{x) = - / f{S)sm{26{L-x))d5., 



(18) 



and the matter density profile inversion problem is solved. One problem is 
very obvious from Eq. (|18|) : One needs to know f{S) in the whole interval 
< i? < oo which is practically impossible. The authors of Ref. Q suggest 
an iteration method to solve this problem. Additional challenges are statis- 
tics and a finite energy resolution, which is "washing out" the edges in the 
profile. One interesting advantage of using solar or supernova neutrinos is 
the sensitivity to asymmetric profiles, i.e., for mass-flavor oscillations there 
is no degeneracy between one profile and the time-inverted version, which 
otherwise (for flavor-flavor oscillations) can only be resolved by suppressed 
three-flavor effects. 



Supernova neutrinos to spy on the Earth's core 

Unlike solar neutrinos, which are limited to energies below 20 MeV, supernova 
neutrinos from a possible galactic supernova explosion have a high-energy 
tail which is closer to the Earth matter resonance energy. This effect is il- 
lustrated in Fig. (right) which compares the energy spectrum between two 
Super-Kamiokande-like detectors with and without Earth matter effects. It 
is obvious from this figure that the difference between the spectra around the 
peak at ~ 20 MeV is tiny, whereas statistically significant deviations can be 
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Fig. 5 Illustration for the tomography of the Earth's core using supernova neutri- 
nos (left) and difference between the event rate spectra (electron antineutrinos) of 
Di and D2 for a Super-Kamiokande-like detectors (right). The difference between 
these spectra corresponds to ~ 35 for the Earth matter effects, i.e., it is highly 
significant. Figures from Ref. I46j |. 

found at larger energies. "'^ Such a scenario could happen if supernova neutri- 
nos were detected by two similar-sized detectors, one on the Earth's surface 
and with the Earth's core in the line of sight (c/.. Fig. [51 left). Note that the 
supernova neutrinos are detected within a very short time frame <^ 24 hr, 
which means that one would actually obtain a "snapshot" of the Earth's 
interior. As it has been demonstrated in Ref. ^4^, for a galactic supernova 
in the distance D — 10 kpc with an energy release of E — 3 10^^ ergs, two 
megaton-size water Cherenkov detectors could measure the density of the 
Earth's core at the per cent level with a number of challenges: First, the 
Earth's mantle density is assumed to be known at the 2% level. Second, the 
solar oscillation parameters have to be known at the 0.2% level. Third, too 
similar supernova fluxes for the different flavors (similar temperatures) and 
deviations from energy equipartition are unfavorable. And fourth, one has 
to have some knowledge on the flavor composition of the flux, possibly from 
detection of different flavors. 



4.3 Neutrino oscillation tomography with neutrino beams 

We now discuss neutrino oscillation tomography with the "man-made" neu- 
trino beams. Neutrino beams are planned or future neutrino sources using 
accelerators, where the neutrino beam is produced by pion/kaon decays (su- 
perbeams, see, e.q., Refs. by muon decays (neutrino factory, see, 

e.g., RefSjj2aJpJa|), or by the decay of unstable nuclei (/3-Beam, see, e.g., 
Refs. [69tll5Hlal36l |). Neutrino beams have, compared to "natural" neutrino 
sources, the advantage that either flux and flavor composition are well-known, 
or a near detector can be used to improve the knowledge on these quantities 

^ For example, at around 34 MeV and 60 MeV deviations between the two curves 
in Fig. 13 (right) can be identified. The difference between these spectra corresponds 
to Ax^ — 35 for the Earth matter effects, i.e., it is highly significant. 
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as well as on the interaction cross sections. There is, however, one major 
obstacle common to all of these experiments: Matter effects especially enter 
in the v^. <-> flavor transition which is suppressed by the small mixing 
angle sin^(20i3). Up to now, this mixing angle is unknown and only an up- 
per bound exists [3| . Experiments within the coming ten years will reveal if 
sin^(20i3) is suitably large for the apphcations discussed here (for a summary, 
see, e.g., Ref. [ssf'). Therefore, the experiment performance has always to be 
evaluated as function of svY?(2Qyi). In this section, we split the discussion 
into conceptual areas linked to tomography with neutrino beams. 

Positional information for a single baseline 

Interesting questions are discussed in Ref. |56| |: Assume we have a beam 
crossing a cavity with a specific density contrast compared to the surrounding 
matter. Then one wants to know 

— How large has the cavity to be to be detected? 

— Can the position of the cavity be measured and if so, how precisely? 

In Ref. a 500 MeV superbeam is assumed with very luminous 200 000 
events in total. The density in the cavity is assumed to be Ig/cm'^ (water), 
the baseline L = 1000km, and sin^(26'i3) — 0.03, where a smaller number 
of events can be compensated by a larger sin^(26'i3). It turns out that the 
cavity has to be longer than about 100 km to be found and its size can be 
measured to about ±50 km. The most important result is that the position 
of the cavity can be reconstructed ±100 km from a single baseline, which 
is very different from X-ray or absorption tomography. However, there is a 
degeneracy in the position between x and L — x which can be only resolved by 
suppressed three-flavor effects. This example demonstrates already one of the 
basic principles of neutrino oscillation tomography: Positional information is 
available already from a single baseline. 

Resolution of structures and edges 

One can learn about the resolution of structures and edges from the numer- 
ical solution of the matter density inversion problem. In Ref. 55] a (sym- 
metrized) REM profile is reconstructed from a single baseline crossing the 
outer core with 14 degrees of freedom using a genetic algorithm. Naturally, 
there are many degenerate profiles close to the la, 2a, and 3a contours, and 
one cannot show a contour in 14-dimensional parameter space. Therefore, we 
show in Fig. El several "typical" representatives close to the la, 2a, and 3a 
contours for a neutrino factory, where the total number of oscillated events 
is for sin^ (26*13) = 0.1 only about a factor of four above currently discussed 
luminosities (see, e.g., Ref. ,33]). From Fig. one can easily read off that 
such an experiment could, in principle, reconstruct the mantle-core-mantle 
structure of the Earth. However, structures smaller than several hundred 
kilometers cannot be resolved. In addition, the mantle-core boundary cannot 
be resolved at a sufficiently high confidence level from a single baseline. Ana- 
lytically, it has been demonstrated in Ref. 55] that structures much smaller 
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Fig. 6 Examples for reconstructed (symmetric) REM profiles from an extremely 
luminous neutrino factory close to the la (upper row), 2a (middle row), and 3a 
(lower row) contours for L = 11 736km (14 d.o.f.). Figure from Ref. |55j |. 



than the oscillation length in matter cannot be resolved - as one would natu- 
rally expect similar to other quantum mechanical phenomena. In conclusion, 
neutrino oscillations in matter are very sensitive towards average densities 
and the arrangement of structure on the length scale of the oscillation length. 
However, neither can edges nor small structures be precisely resolved. 



Density measurement 



Since we know that neutrino oscillations measure more or less the baseline- 
averaged densities = 1/L p{l)dl over long distances plus some sup- 
pressed interference effects, we can use this to discuss possible applications. 
For example, let us assume that we want to perform a simple one-parameter 
measurement of the Earth's inner core density. Because the Earth's mass is 
fixed, we need to correct the average mantle or outer core density for any 
shift of the inner core density. Note, however, that it is the volume-averaged 
density to be corrected, which means that large shifts in the Earth's inner 
core density cause only very small density corrections in the mantle. This 
example illustrates already one potential strength of neutrino oscillation to- 
mography: Since neutrinos from a "vertical" baseline travel similar distances 
in mantle, core, and inner core, there should be no a priori disadvantage for 
the innermost parts of the Earth. In Ref. 1^ a neutrino factory setup from 
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Fig. 7 The measurement of pic (inner core density) as function of the true 
value of sin'^(2&i3) at the la, 2a, and 3a confidence levels (from light to dark 
shaded regions). For the baselines, L = 2 ■ Re combined with a shorter baseline 
L = 3 000 km to reduce correlations is used. The dashed curves correspond to 
fixing the oscillation paramete rs, i .e., to not taking into account correlations and 
degeneracies. Figure from Ref. [63. 

Ref. |33j | with currently anticipated luminosities was chosen to test this hy- 
pothesis for realistic statistics. In order to measure the oscillation parameters, 
the experiment with L = 2 Re was combined with a L = 3 000 km. The pre- 
cision of the measurement can be found in Fig. [7|as function of sin^(20i3). 
One case easily read off that a per cent level measurement is realistic for 
sin^(20i3) > 0.01. Most importantly, the application survives the unknown 
oscillation parameters and the performance is already close to the optimum 
(dashed curves). For smaller values of 0.001 < sin^(26'i3) < 0.01, the corre- 
lations would be much worse without the L = 3 000 km baseline. For large 
values of sin^(2^^i3) > 0.01, the vertical baseline alone is hardly affected by 
correlations with the oscillation parameters: As illustrated in Ref. [B^, CP 
effects are suppressed for very long baselines. Since there is only a number 
of potential high-energy laboratories around the world which could host a 
neutrino factory, we show in Fig. |H1 some examples and the corresponding 
outer and inner core crossing baselines. Obviously, there are potential detec- 
tor locations for some of the laboratories, which are, however, not exactly 
on the L = 2i?£;-axis. Relaxing this baseline constraint somewhat, one can 
show that one can find detector locations for a small drop in precision 67j. 
In summary, this application illustrates that a density measurement could 
be performed with a) reasonable statistics, b) including the correlations with 
the oscillation parameters, and c) reasonably small values of sin^(26'i3). In 
the future, it has to be clarified how large the additional effort for such a 
facility (the vertical storage ring) would be. Note, however, that there are 
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Fig. 8 Positions of three of the major potential neutrino factory laboratories, 
(typical) L — 3 000km detector sites (dashed curves), as well as potential detector 
sites with outer core crossing baselines (below thin solid curves), and inner core 
crossing baselines (within thick solid curves) . The colors of the curves represent the 
different laboratories. Figure from Ref. W^- 



plenty of other applications of a "very long" neutrino factory baseline, such 
as the "magic baseline" to resolve degeneracies (32j {L ~ 7 500 km), the test 
of the MSW effect for sin^(26'i3) = ojel (L > 5 500 km), the mass hierarchy 
measurement for sm^{29i£) = [slt l30l| (L ~ 6 000 km), and the test of the 
"parametric resonance" [aH^ {L » 10 665 km). 

In addition to the described neutrino sources, note that tomography com- 
paring the neutrino and antineutrino disappearance information from atmo- 
spheric neutrinos might, in principle, be possible as well [2^. 



5 Other geophysical aspects of neutrino oscillations 

It is well known that matter density uncertainties spoil the extraction of the 
oscillation parameters from the measurements (see Refs. 0;|27; 62; 22i IsS 
|H3;^2IE3 references therein). In particular for baselines sensitive to Scp, 
such as L ~ 3 000 km at a neutrino factory, the additional correlation with 
the matter density affects the precision measurements of sin^(20i3) and Sep, 
and the CP violation sensitivity. This effect is illustrated in Fig. IHl for the 
precision of ^cp ^-nd different matter density uncertainties Ap* . Especially for 
large sin^(20i3), any uncertainty larger than about 1% affects the precision 
severely. Note that the baseline used for Fig.^is L — 3 000 km, which means 
that the neutrinos travel in an average depth of ^ 120 km up to a maximum 
depth of ^ 180 km. In these depths, the uncertainty among geophysics models 
is currently at the level of 5% [23|. Since the matter density uncertainties may 
affect the competitiveness of a neutrino factory with a superbeam (operated 
at shorter baselines) for large values of sin^(20i3), improved knowledge for 
specifically chosen baselines would be very helpful. 
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Fig. 9 The precision ol the measurement of 5cp for a neutrino factory and the 
simulated value 5cp = 90° as a function of the true value of sin^(2&i3) at the 
la confidence level. The different curves correspond to different allowed matter 
density uncertainties Ap* as described in the plot legend, especially the thick curves 
correspond to no matter density uncertainty (light thick curve) and the often used 
uncertainty Ap* = 5% ■ p (dark thick curve). Figure from Ref. ^5^■ 



6 Summary and conclusions 

In summary, neutrino tomography might be a very complementary approach 
to geophysical methods. For example, neutrinos travel on straight lines with 
almost no uncertainty in their path. Furthermore, neutrino tomography is 
either sensitive to the nucleon density (absorption tomography) or electron 
density (oscillation tomography). In comparison, the paths of seismic waves 
are curved, and there is some uncertainty in them. In addition, the matter 
density has to be reconstructed from the propagation velocity profile by the 
equation of state. This means that neutrino tomography might be a more 
"direct" handle on the matter density and could be very useful to investigate 
specifically localized regions, such as in the lower mantle. Moreover, there is 
no principle reason to prevent neutrinos from penetrating the Earth's core, 
whereas seismic waves are partially reflected at the mantle-core and outer- 
inner core boundaries. Note that though the most precise information on 
deviations from the REM ("Reference Earth Model") in the Earth's mantle 
comes from seismic waves, there are other geophysical methods which might 
be more directly sensitive towards the matter density, such as normal modes, 
mass, and rotational inertia of the Earth. Nevertheless, none of those could 
provide a measurement along a very specific path. 

The main challenges for neutrino tomography might be the existence of 
high-energy neutrino sources for absorption tomography, and the statistics for 
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oscillation tomography. For example, neutrino oscillation tomography could, 
in principle, reconstruct the matter density profile along a single baseline 
due to interference effects among different matter density layers. Note, how- 
ever, that neutrino oscillations are to first order sensitive towards densities 
averaged over the scale of the oscillation length, which means that such so- 
phisticated applications require extremely large luminosities (detector mass 
X source power x running time) and might be very challenging. On the other 
hand, very simple questions, such as a one-parameter measurement of the av- 
erage density along the path or the discrimination between two very specific 
degenerate geophysical models might be feasible within the next decades. 
For example, the achievable precision for the inner core density of the Earth 
with a neutrino factory experiment might be quite comparable (±0.23 g/cm^ 
for sin2(26'i3) = 0.01 and ±0.06 g/cm^ for sin2(26»i3) = 0.1 to current 
precisions given for the density jurnp at the inner-core boundary from geo- 
physics {e.g., ±0.18 g/cm^ in Ref. |43). We therefore conclude that it will 
be important that the right and simple questions be asked by discussions 
between neutrino physicists and geophysicists. 
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